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SECTION  1 
INTRODUCTION 


Ground  wave  propagation  is  possible  at  all  radio  frequencies  but 
often  neglected  in  the  HF  band  (3-30  MHz)  because  of  the  much  longer  paths 
possible  with  sky  wave  propagation.  Typical  ground  wave  propagation  paths 
are  measured  in  tens  to  hundreds  of  kilometers  while  typical  sky  wave 
paths  extend  to  thousands  of  kilometers.  The  purpose  of  this  paper  is  to 
consider  various  aspects  of  HF  ground  wave  propagation  over  smooth  and 
irregular  terrain.  Although  HF  frequencies  are  emphasized,  many  of  the 
results  are  applicable  to  a  much  wider  range  of  frequencies. 

A  military  situation  involving  a  large  nuclear  exchange  is  one 
example  where  ground  waves  may  be  considerably  more  effective  than  sky 
waves.  Nuclear  explosions  cause  severe  D-region  absorption  which  can 
reduce  or  eliminate  conventional  HF  sky  wave  modes.  These  explosions  are 
expected  to  have  little  or  no  effect  upon  ground  waves.  While  sky  waves 
can  be  degraded  by  attacking  the  propagation  medium  (the  ionosphere)  such 
a  tactic  is  nearly  impossible  with  ground  waves.  Ground  waves  propagate 
along  the  ground/air  interface  and  are  only  sensitive  to  substantial 
changes  in  this  interface.  Even  a  large  nuclear  exchange  would  not  be 
expected  to  greatly  alter  the  basic  topography  of  a  region  and  hence,  the 
ground  wave  propagation. 

On  the  other  hand,  ground  wave  propagation  is  vulnerable  to  jam¬ 
ming  or  to  direct  destruction  of  the  terminals  as  is  sky  wave  propagation. 
However,  jamming  of  ground  wave  circuits  can  be  made  difficult  by  taking 
advantage  of  the  limited  range  of  ground  waves.  If  sky  wave  propagation 
is  suppressed  by  proper  choice  of  frequencies  and  antennas,  then  an  enemy 
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must  place  a  jammer  near  or  within  sight  of  the  terminals  being  used  to 
effectively  jam  them. 

Ground  wave  propagation  generally  requires  as  much  or  more  trans¬ 
mitter  power  and  as  good  or  better  antennas  than  sky  wave  propagation. 
Higher  transmitter  power  enhances  enemy  direction  finding  and  eavesdrop¬ 
ping  but  is  compensated  by  the  relatively  small  ground  wave  coverage 
regions.  Because  ground  wave  propagation  is  better  at  lower  frequencies, 
it  is  tempting  to  consider  operation  at  the  low  end  of  the  HF  band  or  even 
into  the  MF.  However,  efficient  antennas  at  these  frequencies  are  large 
and  probably  impractical  for  certain  mobile  applications.  For  these 
applications,  it  may  be  better  to  operate  at  middle  HF  frequencies  where 
efficient  antennas  of  moderate  size  are  possible. 

In  general,  ground  wave  propagation  offers  worthwhile  capabili¬ 
ties  for  short  distance  military  communications  systems  intended  to  sur¬ 
vive  a  nuclear  attack.  Network  communications  systems  with  many  closely 
spaced  nodes  may  find  ground  waves  attractive. 

Because  the  subject  of  ground  wave  propagation  is  broad  and 
complex,  this  report  concentrates  on  those  aspects  which  have  recently 
been  studied  for  DNA.  An  extensive  survey  of  the  ground  wave  literature 
has  been  made  leading  to  the  selection  of  a  theoretical  formalism  appro¬ 
priate  to  computer  coding.  The  formalism  chosen  was  developed  by  Hendri- 
fcus  Bremmer  over  30  years  ago  and  applies  to  propagation  over  a  smooth 
spherical  earth.  This  formalism  has  been  coded  in  Fortran  IV  and  is 
reported  here.  The  computer  routine  has  been  kept  relatively  simple  and 
does  not  yet  include  corrections  for  atmospheric  refraction  or  for  rough 
earth  cases.  These  subjects  are,  however,  discussed  in  some  detail. 
Finally,  the  lack  of  experimental  data  is  discussed  along  with  a  sugges¬ 
tion  for  a  future  experiment. 

This  work  was  supported  by  two  DNA  contracts:  DNA  001-80-C-0022 
and  DNA  001-80-C-0225. 


6 


SECTION  2 

GENERAL  CONSIDERATIONS 


There  is  often  confusion  as  to  the  physical  reality  of  ground 
waves  (or  surface  or  terrestrial  waves).  They  are  that  portion  of  elec¬ 
tromagnetic  waves  propagated  through  space  and  affected  by  the  presence  of 
the  ground.  They  do  not  include  any  portion  of  the  wave  reflected  from 
anything  other  than  the  ground.  For  example,  ionospheric  sky  waves  or 
tropospheric  waves  are  separate  phenomena.  Line-of-sight  propagation 
between  two  spacecraft  in  empty  space  away  from  the  earth  is  also  a  separ¬ 
ate  phenomenon  although  "line-of-sight"  propagation  between  two  aircraft 
aloft  is  technically  ground  wave  propagation.  As  the  aircraft  fly  higher, 
ground  wave  propagation  approaches  line-of-sight  propagation.  In  other 
words,  the  dielectric  properties  (dielectric  constant  and  conductivity)  as 
well  as  the  topology  of  the  ground  influence  the  propagation  of  electro¬ 
magnetic  waves  "near"  the  ground.  "Near"  is  defined  in  terms  of  the 
number  of  wavelengths  that  the  transmitter  and  receiver  are  above  the 
ground.  At  two  meters  above  the  ground,  visible  light  is  hardly  affected 
at  all  by  the  ground  while  HF  transmissions  are  strongly  affected. 

Ground  wave  propagation  is  a  diffraction  effect  but  not  one  of 
the  simpler  diffraction  effects  often  considered  in  the  study  of  optics. 2 
Although  the  general  concept  of  diffraction  around  an  edge  applies  to 
ground  wave  propagation  around  the  earth,  the  mathematical  formalism  is 
considerably  more  complex  for  ground  wave  propagation.  At  one  time,  this 
subject  was  considered  to  be  among  the  most  difficult  in  theoretical 
physics.  To  obtain  any  solution  at  all  requires  several  assumptions. 
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Arnold  Sommerfeld  solved  the  general  problem  of  radio  propaga¬ 
tion  from  a  short  vertical  antenna  over  a  finitely  conducting  plane  earth 
in  1909.3  The  rigorous  results  obtained  by  Sommerfeld  were  applicable 
to  short  transmission  distances  where  the  curvature  of  the  earth  could  be 
ignored.  It  was  already  known  that  the  solution  for  a  spherical  earth 
could  be  obtained  in  terms  of  a  series  of  zonal  harmonics.  Unfortunately, 
this  series  converged  so  poorly  that  no  practical  results  could  be  obtain¬ 
ed  for  radio  waves. 

In  1918,  Watson^  developed  a  transformation  with  the  aid  of  an 
integral  in  the  complex  plane  which  transformed  the  solution  in  terms  of 
zonal  harmonics  into  a  residue  series  which  converged  rapidly  and  could  be 
evaluated  numerically.  From  1926  through  1931,  Sormierfeld,  Van  der  Pol, 
Niessen,  and  Wise  obtained  several  independent  solutions  of  the  problem 
using  Watson's  transformation.  During  this  period,  a  fundamental  error 
was  also  uncovered  and  corrected  in  Sommerfeld's  original  work.  In  the 
mid-1930's,  Bremmer  and  Van  der  Pol  in  Europe  and  Norton  in  the  United 
States  extended  the  early  work  to  the  point  where  it  became  practical  to 
use.5»6,7,8,9,10  They  simplified  the  mathematics,  performed  extensive 
numerical  computations,  and  helped  to  better  define  the  physical  nature  of 
ground  waves.  After  the  Second  World  War,  many  other  investigators 
considered  various  aspects  of  the  problem  such  as  an  inhomogeneous  or 
rough  earth.  The  most  generally  useful  results,  however,  are  still  those 
of  Bremmer  or  Norton.  We  have  arbitrarily  chosen  to  follow  Bremmer' s 
formal  ism*  because  it  appears  to  be  the  easiest  to  use. 

Bremmer's  calculations  reveal  several  general  characteristics  of 
ground  waves.  For  instance,  they  propagate  much  better  over  water  than 
over  land,  and  better  over  wet  land  than  dry  land  because  the  high  conduc¬ 
tivity  of  water  and  wet  land  produces  less  absorption.  Ground  waves  are 
more  influenced  by  the  conductivity  of  the  ground  (soil  type)  than  by  the 
curvature  of  the  earth  (for  transmitter  and  receiver  on  the  earth).  They 
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are  not  abruptly  changed  by  the  horizon.  As  with  most  diffraction 
processes,  there  are  no  distinct  shadows  behind  obstacles  such  as  the 
horizon.  Also,  long  wavelengths  propagate  farther  around  the  earth  than 
do  short  wavelengths.  Ground  waves  are  usually  vertically  polarized 
because  horizontal  polarization  is  readily  short-circuited  by  the  earth. 
Hence,  horizontally  polarized  antennas  such  as  might  be  used  for  sky  wave 
propagation  are  relatively  ineffective  for  ground  wave  propagation. 
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SECTION  3 

BREMMER'S  EQUATIONS/GROUND  WAVE  PROGRAM 


This  section  presents  a  summary  of  Bremmer's  approach  to  ground 
wave  theory  and  shows  how  this  theory  has  been  used  to  create  a  simple 
computer  routine  for  computing  ground  wave  field  intensities  around  a 
smooth  spherical  earth. 


As  mentioned  in  the  last  section,  the  rigorous  theory  for  wave 
propagation  around  a  sphere  is  a  special  application  of  diffraction 
theory.  It  is  mathematically  difficult  because  of  the  finite  size  of  the 
sphere.  Much  simpler  solutions  are  possible  for  a  very  small  or  large 
sphere.  The  general  solution  in  terms  of  a  series  of  zonal  harmonics  is 
given  by^ 


H  = 


IL 

cb 


.V 


*  "<o  I <2ntl>Rn 
n=0 


xry 


(koa) 


41}(kob)  pn{cose)] 


(i) 


where  H  is  the  magnitude  of  Hertzian  vector.  This  Hertzian  vector  is 
related  to  the  conventional  vector  potential.  A,  by  the  relation 


n 


A  =  1  3A 


(2) 


From  a  knowledge  of  the  Hertzian  vector,  the  electric  and  magnetic  field 
components  follow  directly. 
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The  effective  reflection  coefficient  in  (1)  is  defined  by 


1ST  "  1”n<X)1^1a 
-I|.»  [xv.)]  x=kia 


and  the  functions  cn( :)(x)  and  Tn(x)  are  given  by 


»  (x)  =  JT  3  <x)  ,  x"  (-i  U  "  (llii) 

n  V  2x  n+1/2  ^  x  dx ;  v  x 


where  H  and  J  are  the  conventional  Hankel  and  Bessel  functions.  Other 
quantities  in  (1)  are  defined  as 

IL  =  electric  dipole  moment  of  the  assumed  small 
vertical  dipole  antenna 

s  =  distance  between  transmitter  and  receiver 

b  =  radial  distance  of  antenna  from  center  of  assumed 
spherical  coordinate  system 

k0  , k x  =  radial  wave  numbers  outside  and  inside  the  earth 
a  =  radius  of  sphere  (in  this  case,  the  earth). 


The  numerical  evaluation  of  (1)  is  very  difficult  and  for  a  long 
time  remained  unknown.  The  problem  is  the  convergence  of  the  series  for  a 
finite  radius,  a.  If  a  is  small  compared  to  the  wavelength,  the  series 
converges  rapidly,  and  the  first  term  yields  Rayleigh  scattering.  The 
greater  the  value  of  the  parameter  kxa  =  2*a/x,  the  more  slowly  the  series 
converges.  The  most  significant  terms  are  those  of  order  n  -  kja.  In  the 
radio  case,  kLa  varies  from  103  to  109  requiring  the  summation  of  a  great 
number  of  terms. 


Because  such  a  summation  is  impractical,  Watson  developed  a 
transformation  which  converted  (1)  into  a  contour  integral  in  the  complex 
n-plane.  This  reduced  the  problem  to  a  summation  of  residues  which  can  be 
mastered  numerically.  Bremmer  introduced  approximations  for  the  resulting 
Hankel  functions  and  zonal  harmonics  leading  to  a  simpler  expression  for 
the  r.m.s.  value  of  the  electric  field  at  a  distance,  D^,  from  a  short 
vertical  dipole  antenna. 1  That  expression  is  shown  below  with  two 
significant  modifications  to  bring  these  computations  in  line  with  the 
conventions  of  the  CCIR  (see  Appendix).  The  electric  field  has  been 
normalized  to  give  a  value  of  346.4  mv/m  at  1  km  if  the  assumed  dipole 
transmitting  antenna  is  placed  on  a  perfectly  conducting  plane  earth.  The 
same  dipole  antenna  will  give  a  field  of  173.2  mv/m  at  1  km  in  free  space 
which  is  the  free  space  field  from  a  1  kw  isotropic  antenna.  Bremmer  and 
others  use  300  mv/m  instead  of  346.4  which  corresponds  to  the  field  from  a 
dipole  antenna  radiating  1  kw  over  a  perfectly  conducting  plane  earth. 
The  antenna  assumed  here  radiates  4/3  kw  under  these  conditions.  Another 
modification  to  Bremmer1 s  equations  involves  the  distance  Dkm .  While 
Bremmer  uses  straight  line  and  great  circle  distance  equivalently,  the 
following  expressions  use  the  latter  only. 


E  =  [ 346 -4-y  2-itx  /"V  (h1+a)z  +  (h2+a)z  -  2(h1+a)  (h2+a)  cos  (D^/a ) ] 


x  I  I  MM  MM  e1  TsX/(2t_-1/62)  I  mv/m  (6) 

s=0  5,1 

where 

D.  =  distance  along  the  surface  of  the  earth  from  transmit¬ 
ter  to  receiver  in  km. 


x  •  -053693  0^/^1/3  (7) 

xm  =  wavelength  in  m 
f s ( h ! )  =  height  gain  factor 
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(8) 


x^-2,s\l/2 

.  / 


Hl/3(l)  [i  (*,2-2ts)3/2) 

H,/3(,)  (i  (-2ts)3/2] 


Xl2  =  .03674  h,/x2/3 
*1  1  m 

hx  =  height  of  transmitter  above  ground  in  m 

h2  =  height  of  receiver  in  m 

6  ,  K  i(135'-ifpe) 

e  e 


.002924  X, 


1/3  >/e2  +  3-6xl°25qe 
tyle- 1)  2  +  3. 6x10 2 5o2  X2 


*e  =  tan 


-l 


E  \  1  .  _  1  f  £“  1 

- )  —  tan  ( - 

6x10 12  2  6x10 12  c 


e  m 


\n 


The  values  of  xs  follow  from 

t  =  Re  t  +  i  Im  t 
s  s  s 

and  assuming  Ke  is  small: 

Re  t0  =  .928  +  K  cos(45*  +  *  )  +  1.237  K  3  cos(75*  +  3*  ) 

y  6  6  G  6 

-  .5  K  k  cos (4  ¥  )  -  2.755  K  5  cos(75°  -  5v  )  ... 
e  e  e  e 

Re  TX  =  1.622  +  Kgcos(45°  +  *  )  +  2.163  Kg3  cos(75°  +  3ffi) 

-  .5  Ke4  cos(4Te)  -  8.422  Kg5  cos(75*  -  5Ye)  ... 


(9) 

(10) 


(ID 


(12) 


13 


Re  t2  =  2.191  +  Kg  cos(45°  +  +  2.921  Kg3  cos  (75"  +  3yg) 

-  .5  Kg4  cos(4*e)  -  15.36  Kg5  cos(75°  -  5yg)  ... 

Re  t3  =  2.694  +  K  cos  (45°  +  y  )  +  3.592  K  3  cos  (75*  +  3y  ) 

-  .5  Kg4  cos  (4»e)  -  23.227  Kg5  cos  (75°  -  5 f  )  ... 

Re  t  =  1.116  (s  +  3/4) 2/ 3  +  K  cos  (45°  +  y  ) 

b  e  6 

+  1.488  (s  +  3/4)  2/3  Kg3  cos  (75°  +  3yg) 

-  .5  Kg4  cos  (4yg) 

-  3.987  (s  +  3/4) *♦/ 3  Kg5  cos  (75°  -  5yg)  ...  for  s  >  3 

Im  tq  =  1.607  -  Ke  sin  (45°  +  yg)  -  1.237  Kg3  sin  (75°  +  3fg) 

+  .5  K  4  cos  (4?)  -  2.755  K  5  cos  (75°  -  5*J  ... 

Im  Tl  =  2.810  -  Kg  sin  (45°  +  yg)  -  2.163  Kg3  sin  (75°  +  3yg) 

+  .5  K  4  sin  (4y  )  -  8.422  K  5  sin  (75°  -  5 ... 

Im  t2  =  3.795  -  K  sin  (45°  +  »  )  -  2.921  K  3  cos  (75°  +  3y  ) 

+  .5  Kg4  sin  (4*e)  -  15.36  Kg5  sin  (75°  -  5^)  ... 

Im  t 3  =  4.663  +  K  sin  (45°  +  ?  )  +  3.592  K  3  cos  (75°  +  3y  ) 

“  6  6  0 

+  .5  Ke4  sin  (4yg)  -  23.227  Kg5  sin  (75°  -  5yg)  ... 
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Im 


1  - 932 ( s  +  3/4) 2/ 3  .  K  sin  (45°  +  T  ) 


1 . 488 ( s  +  3/4) 2/3  k  3  sin  (750  +  3y  \ 
e  e 


+  .5  K  4  sin  (Ay  ) 
e  e 

-  3.987(s  +  yn)"/*  Kg5  sin  (75°  -  5»  )  ... 


for  s  >  3 


For  Ke  large: 


n  cos  (15°  -  »e)  sin(2vp) 

Re  t  o  =  .4043  +  .6183  _ !_  -  .2364  —  —  e 


K, 


-  .0533 


cos  (15°  +  3¥e) 


Ke3 


+  .00226 


Ke2 

cos  (60°  -  4Ye) 


Re  T1  -  1.288  +  .194  1°.!-  (158  '  ^  -  .0073 

Ke  Ke2 


cos  (15°  -  3vp ) 


---  -ip /  cos  (60°  -  4vp) 

+  .0120  - £_  -  .00160  - - - 


K, 


Ke4 


Re 


T  _  * 


1.116  (s  +  1/4) 2/3  +  _ -2241  C°S  (15  I  4'e) 

(S  +  174)2/3  Ke 


Im 


To 


for  s  >  1 


sin  (15°  -  <Fe)  cos  ZVp 

- —  +  .2364  - ® 


,7003  -  .6183  - : - —  +  .2364 

Kq 


Ke2 


-  .0533 


sin  (15°  +  3ip) 


Ke3 


e/  sin  (60°  -  40 

—  -  .00226  -  e' 


Kp4 


15 


2.232  -  .1940 


+  .0073 


Im  ti 


Im  Ts 


sin  (15°  -  fj 


sin  (15°  +  3fP) 

+  .0120  _ —  +  . 


xe3 


00160 


1.932  (s  +  1/4) 


2/  3 


.2241 


(s  +  1/4) 2/3 


cos  (2^) 

"  Ke2 
sin  (60*  -  4»e) 

*7 

sin  (15*  -  fe) 

■  ■  -  ■'  -»■■■-  •  •  • 

Ke  for  s>l 


These  expressions  apply  to  a  vertical  dipole  transmitting 
antenna  and  to  a  separation  between  transmitter  and  receiver  of  at  least 

Dkm  ^  ^  Am  1^3* 

For  a  horizontal  dipole  antenna,  the  above  equations  can  be  used 
(in  the  direction  of  the  maximum  field)  providing  that  6e,  Ke,  and  i|>e  are 
replaced  by  sm,  and  90’  -  ¥m.  These  latter  quantities  are  defined  as 
follows: 


=  K,  ei(45#  + 

(13) 

\  i/3 

=  .002924  ^ 

—  (14) 

*W(e  -  l)  2  +  3.6  x  1025  ae2  X,,,2 


1  tan-1  ( _ _ ) 

2  6  x  10lz  ^ 


(15) 


It  should  also  be  noted  that  ae  is  defined  in  terms  of  e.m.u. 
following  Bremmer's  preference.  In  terms  of  more  customary  usage,  this 
can  be  written  as 

oQ  =  1011  oo 

e  a 

e.m.u.  mhos/m 
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The  height  gain  factors  in  (6)  and  (8)  are  made  up  of  Hankel 
functions.  Using  Bremmer's  approximation  for  these  functions,  they  can  be 
written  as 


a)  for  hlm  >  50  (approximately): 

fs(M  -  je1'  i""'  *  T  'S'  '  2^)3/2'  „  “6) 

5e  V  -  ( 

,  .2038  .3342 

(x,2-  2ts)2/2'  (xt2-  2,$)5  ‘ 


where 


.[in/4  -  ~  (x  2  -  2t  )3/2] 


1  +  i 


.2083 


k,‘ 


2  )3/2  ■ 


A0  =  .3582  ei  120' 

Ay  =  .3129  e-i  60* 

A2  =  .2903  e1  120' 

A3  =  .2760  e"1  60# 

A  =  .3440  - e'iir/3 

S  (s  +  3/4) 


(17) 


for  s  >  3 


b)  for  h  <  50  Am2^3  (approximately) : 


Vhi > 


1  +  6.283  ( - ! 

xi/3 


39.48 


(l-x2/3  6e  ts) 
x-/3  6e 


(18) 


where 


x 


4xl07 


X 


m 
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The  height  gain  factor  for  the  receiver  fs(h2)  can  be  computed 
as  above  by  substituting  h2  for  hr.  Similarly,  the  height  gain  factors 
for  a  horizontal  dipole  can  be  computed  from  the  equations  above  by 
substituting  ^  for  6e.  Bremmer  maintains  that  methods  (a)  and  (b) 
above  are  sometimes  equally  suitable. 

Using  equations  (6)  through  (18),  a  ground  wave  program  has  been 
constructed  and  is  given  in  Table  1.  By  modifying  the  input  and  output, 
this  program  has  also  been  converted  into  a  subroutine  for  the  nuclear 
effects  code  known  as  HFNET.  The  necessary  inputs  (including  units)  are 
shown  in  the  program  listing.  This  program  applies  to  ground  wave  propa¬ 
gation  around  a  smooth  homogeneous  spherical  earth.  It  does  not  contain 
any  provisions  for  a  rough  or  inhomogeneous  earth.  It  also  does  not 
contain  the  effects  of  atmospheric  refraction.  Despite  these  limitations, 
the  program  accurately  computes  ground  wave  field  strengths  in  many  prac¬ 
tical  cases  where  there  are  no  large  departures  from  the  assumptions  built 
into  the  program.  The  program  is  not  limited  to  the  HF  band  but  applies 
to  a  wide  range  of  radio  frequencies.  It  has  been  checked  out  from  15  KHz 
to  100  GHz  but  may  be  applicable  to  an  even  wider  range  of  frequencies. 
It  is  also  applicable  to  a  wide  range  of  separations  between  transmitter 
and  receiver.  As  will  be  shown  below,  the  routine  does  fail  under  some 
circumstances  where  the  transmitter  and  receiver  are  close  together  (in 
terms  of  number  of  wavelengths)  or  high  above  the  earth  in  full  view  of 
each  other.  Under  these  conditions,  a  "geometric  optical"  model  is  more 
appropriate  and  will  be  added  to  the  program  in  the  future.  Although  the 
residue  series  model  used  in  this  program  is  relatively  less  efficient 
when  the  receiver  is  in  view  of  the  transmitter,  the  results  are  still 
accurate  except  as  noted.  Many  more  terms  in  the  residue  series  need  to 
be  summed  *rf>en  the  receiver  is  above  the  transmitter's  geometrical  hori¬ 
zon. 

Extensive  checkout  of  the  ground  wave  program  has  been  accom¬ 
plished  and  is  shown  in  Figures  1  through  8.  (All  of  this  checkout  was 


000 1 
0002 
0003 
0004 

ooos 

0006 

0007 

0008 

0009 

0010 

0011 

0017 

0013 


0014 


0015 


0016 


0017 
0018 
0019 
0020 
002 1 
0027 
0023 
0024 

0025 

0026 

0027 

0028 

0029 

0030 

003! 

0037 

0033 


0034 


Table  1.  Ground  wave  program  listing. 

C  PROGRAM i  GROUND  WAVE 

C 

C  **** ************************************ 

c 

C  BREMMER  GROUND  WAVE  MODEL 

C 

C  CODED  BY  GORDON  J.  FOLKS  4/80 

C 

c  **************************************** 

DIMENSION  TAUR< 0 J 3) r  TAUI<0:3>*  HT<2>*  CSQ<2> 

REAL  K»K2fK3rK4*K5r LAMBDA 

COMPLEX  DELTA*  J*  A<013>,  DIF*  FI*  F2*  AX 
COMPLEX  SUM*  F3*  SUMX*  SSDIF*  HGF<  2) *  TAU 
OPEN  <UNIT“1 * NAhE= 'GND. DAT ' * TYPE* 'OLD' *  READONLY  > 

OPEN  ( UNI  I =2*  NAME” ' GND • OUT ' *  TYPE* ' NEW ' ) 

Cl  *  1./3. 

C2  -  2.  *  Cl 
C3  *  4 .  *  Cl 
C4  -  Cl  /  2. 

J  «  CMPLX ( 0 . *  1  ♦  ) 

PI  *  3.1415976535 
RAD  a  6.371E3 

C  **************************************** 

c 

r  INPUT  SECTION 

C 

c  **************************************** 

C 

C  READ  TRANSMITTER  POWER  IN  KILOWATTS*  TRANSMITTER  WAVELENGTH 

C  IN  METERS*  AND  GREAT  CIRCLE  DISTANCE  ALONG  THE  SURFACE  OF 

l-  THE  EARTH  FROM  TRANSMITTER  TO  RECEIVER  IN  KIl  OMETERS 

C 

READ  (1**)  PWR*  LAMBDA,  D 

C 

t:  READ  HEIGHTS  OF  TRANSMITTER  AND  RECEIVER  ABOVE  THE  GROUND 

C  IN  METERS  AND  WHETHER  THEY  ARE  VERTICALLY  OR  HORIZONTALLY 

C  POLAR I 7ED  <0  -  VERTICAL*  1  -  HORIZONTAL) 

C 

READ  <1**>  HT*  NPOL 

C 

C  READ  GROUND  CHARACTERISTICS 

C  EPSLN  -  RELATIVE  DIELECTRIC  CONSTANT 

C  SIGMA  a  CONDUCTIVITY  IN  MHOS/METER 

C 

REAP  (1**)  EPSLN*  STGMA 

C  *************************************** 

c 

C  COMPUTATION  SECTION 

c 

c  *************************************** 

SIGMA  =  SIGMA  *  I.E-U 

CHI  a  .053693  *  D  /  L AMBDA**C1 

HI  R  HT  < 1 )  *  t.E-3 

H2  =  HT < 2 )  *  l.E-3 

DIS  a  (HI  +  RAD ) **2  +  ( H2  I  RAD)*#2 

DIS  a  SORT <  DIS  -  2*<H1  +  RAD)  *  <M2  f  RAD)  *  COS <  D/RAD ) ) 

K  =  .002924  *  LAMRDA**Cl 

K  «  K  /  SQR I ( SORT (< EPSLN  -  1.)**?  +  3.6E25  *  SIGMA**2 
1  *  LAMBDA**?)) 

TF  (NPOL.FO.t)  GO  TO  50 

K  *  K  *  SQR  I  <  EF'SL  N#*2  f  3.6E25  *  SIGMA**2  *  L  AMBDA**2 ) 

FSI  ^  A TAN (EPS!  N/(6. E12  *  SIGMA  *  LAMBDA)) 

PSI  *  FSI  -  .Ti  »  ATAN< (EPSLN  -  I.)  /  (6.E12  *  SIGMA  *  LAMBDA)) 
DELIA  -  K  4  FXE(.I  *  ((3.  *  El  /  4.)  -  ESI)) 

GO  in  100 

50  ESI  «  .5  *  AIAN< (EPSLN  -  1.)  /  (A.E12  *  SIGMA  *  LAMBDA)) 

DELTA  =  K  »  EXE ( J  *  (<PI  /  4.)  +  ESI)) 

ESI  *  El  /  2.  -  PSI 
C 

E  RESIDUE  SERIFS 

C 

100  K2  ^  K»*2 
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Table  1.  (Cont.) 


0035 

K3  K**3 

OO  3  A 

K4  ^  K**4 

003/ 

K5  -  K**5 

003H 

If  <K .  GT .  0  .  A )  GO  TO  200 

00  T9 

ANSI  -  TPI  /  4.)  4-  PS I 

0040 

ANG?  -  75.  *  PI  /  ISO.  +3.  *  PS I 

004  f 

ANG  3  -  4.  *  PS I 

0041’ 

ANG4  =*  75.  *  PT  /  1B0.  -  5.  *  PSI 

0043 

CANG1  -*  COS T ANG 1) 

0044 

SANG 1  -  S1NTANG1) 

0045 

GANG?  *  COST ANG2 ) 

00 4  A 

SANG?  -  S INTANG?) 

004  7 

GANG 3  *  COS  <  ANG3 ) 

0048 

GANG.!  =-  SINTANG3) 

0049 

GANG 4  -  COST ANG4 ) 

0050 

SANG4  -  S 1 N ( ANG4 ) 

0051 

TAURTO)  “  .  92B  4  MtCANGl  4  1 . 237*K3*CANG2  .5*N4*CANG3  - 

l  2.755*K5*CANG4 

0052 

r AUR (  1 )  J.A22  4  N*CANG1  +  2.  163*h3*CANG2  -  .5*N4*CANG3  - 

1  8,422*K5*0ANG4 

005  T 

TAUR( 2 )  =  2.191  4  KfCANGl  4  2 . 921 *K3*C ANG2  -  .5*K4*CANG3  - 
1  15.3A*K5*CANG4 

0054 

T AUR ( 3 )  “  ? , A 9 3 7  4  K*CANG1  4  3 . 5919*K3*CANG7  -  . 5*K4*0ANG3 

1  73.2?AB*K5*CANG4 

0055 

AUTTO)  -  1.A07  -  K*SANG1  -  1 . 237*K3*SANO?  4  ,5*K4*SANG3  - 
1  2. 755*K5*SANG4 

005/* 

TAHITI)  “  ?.  BIO  -  K*SANG1  -  2 . 1 A3*K3*SANG?  4  .5*N4*SANG3  - 
1  B.422*K5*SAN04 

005  7 

TAUT  T  2)  »  3.795  -  K*SAN0.1  -  2 . 921*K3*SANG?  +  .5*N4*SANG3 
t  15.3A*K5*SANG4 

005ft 

I  AO  1(3)  -  4 . A  A  3  3  -  N*SANG 1  -  3.59t9*K3*SANG?  4  .5*K4*SANG3 
1  23. ?2AB*K5*SANG4 

0059 

GO  10  300 

00  *0 

200 

ANSI  -  15.  *  PT  /  1  BO  *  -  PSI 

OOA1 

ANG?  2.  *  PSI 

00  A? 

ANG3  “  15.  *  PI  / 1B0  4  3.  *  PSI 

oo  A. 5 

ANG 4  -  PI  /  3.  4.  *  PSI 

006  4 

RCANfil  =  COSTANGJ. ) 

00A5 

HSANG1  •»  SIN  T  ANG  1) 

00  A  A 

CANG?  ~  COS TANG?) 

00  A  7 

SANG?  “  SINT ANG?) 

00  Aft 

CANG3  ■=  COST  ANG3  ) 

00A9 

SANG 3  •-*  SI  NT  ANG 3) 

00  70 

CANG4  COS  <  ANG 4  ) 

00  7 1 

SANG 4  -  S  INT  ANG 4 ) 

0072 

fAOR(O)  ==  .4043  4  .  A I  B*HCANG  l /h  *  .  23A4*GANG2/h? 

1  .0533*i:ANL»37K3  4  . 0022A*CANG4/K4 

0073 

I  AUR (  t  )  ="  1.2BB  ♦  .1V4#BCANG1/K  . 00734SANG2/K2  4 

1  . 01 ?0*CANG3/K3  -  . 001 60*CANG4 7K4 

00/4 

1  All  1  <  0  >  -  .  7003  .  A1H3*HSANG1 7K  4  .  23A4*C ANG2/K2  - 

1  .<>533*SANGJ/K3  -  . 00226*SANG4/K 4 

0075 

1  All  1  <  1 )  =  2.237  .  19404BSANG1/K  f  . 0073#CANG2/K?  4 

1  . 0 1 ?0*SANG3/K3  4  . 001 60*SANG47K4 

00  fh 

Ml  730  I  a  2>  3 

00  77 

lAORTl)  «  1.ilA*((  4  . 25>**C2  4  . 224 1 *HCANG1  / 

1  TK*<I  4  .  25  )  *  #(;? ) 

00  /B 

TAO  1(1)  1  » 9 3 2 * <  I  4  .?5)**C?  -  .  224  1  HBSANGl  / 

1  <K*TI  4  . 25 )**C2) 

00  7  V 

:mo 

CONI INOb 

C 

c 

300 

HUGH!  GAIN  FACTORS 

00  BO 

COM  --  50.  *  l  AHttUA**C2 

00  Bt 

sim  =-•  o. 

0002 

E  -  BAH. 3  *  SORT (CM I)  *  1.K3/HIS 

0003 

ATO)  “  .3582  *  bXPTJ  *  C2  *  PI) 

00B4 

ATI)  *  .3179  *  EXP!  .1  *  Cl  *PI) 

00B5 

AT  2)  -  .2903  *  bXPTJ  *  C2  *  PI) 

OOBA 

AT  3)  =  .2  760  *  EXPT-.J  *  Cl  ¥PT> 

00  B  7 

Ml  399  N  -  1 » 2 

008B 

('SOT N )  =  . 03 A 74  *  Ml  < N  >  7  1  AMhI«A**C2 

OOBV 

3V9 

(ONI INOb 
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0090 
0091 
0092 
00 93 
0094 
0095 
0096 

009/ 

0098 

0099 

0100 

0101 
0102 
0103 
0104 
0105 
0106 
010/ 
0108 
0109 
0110 
0111 
0112 
0113 
01  14 
0115 
01  16 
01  17 
Ot  18 

01  19 

0 1  20 
0121 

0122 

0123 

0124 

0125 

0126 


OJ  2/ 
0128 
0129 
0130 


0131 

0132 

0133 


0134 


410 


450 


460 

500 


510 

530 


550 


/oo 


899 

900 

r: 

c 

c 

c 

r 


r: 

c 

c 


K*CANG1  f  1 .4881*<H  .75)**C2 
3 . 986/* ( 1+ . 75 ) **C3*K5*CANG4 
K*SANG1  -  1 • 4881* ( If.75)*#C2 
3* 9867* ( If. 75 > **C3*K5*SANG4 


HO  899  I  =  0 » 500 
If-  ( 1  •  G I  .  3 )  GO  TO  410 
AX  “  A ( 1 ) 

TAO  =  CMPLX  (  TAllf\'(  I )  » TAIII  ( I  )  ) 

GO  TO  500 

If  (K.GT.0.6)  GO  TO  450 
TAllKX  =  1 , 116*< If ,75)**C2  f 
1  *K3*CANG2  -  ♦ 5*K4*CANG3  - 
1AUIX  -  1 .932*< If .75>**C2  - 
1  *K3*SANG2  f  « 5*K4*SANG3  - 
GO  TO  460 

IAIJRX  =  1  . 11A*<  If  .2b)**C2  f  .  224 1 *BCANG1  / 

1  (K*< If »25>**C2) 

TAU I X  =  1 . 932*  (If  .25  >**02  -*  .2241*BSANG1  / 

1  <  K*  < If. 25  >  **C2 ) 

AX  *  .3440  *  <~l)**<Ifl>  *  EXP ( - J  *  PI/3 )  /  < If . 75 ) **C4 
TAU  =  CMPLX <  TAURX  t  1AUIX> 

HO  700  N=  1 1 2 

IF  <HT<N) .LT.COM)  GO  TO  550 
HIF  -  CSO<N>  -  2.  *  TAU 

FI  =  -  .1  *  PI/4,  f  .J  *  (  (SORT (DIF)  >**3>  /  3. 

P-  2  -  1 .  -  .1  *  5./ <24.  *  ( SORT  <DIF>)**3)  -  385*/<1152 
F  3  =  1 .  fJ*  5 • / ( 24 .  *  < SORT < DIF > )#*3> 

IF  < ADS < REAL < F 1 ) ) . GT . 89 . >  GO  TO  510 
HUK<N>  =  F2  *  EXF‘< F 1  >  -  F3  *  EXP(-Fl) 

SSDJF  *  SURMSORT  (DIF)  > 

HliMN)  =  MGF(N)  *  AX  /  (Uhl  I A  *  SSD1F) 

GO  10  530 
HOF  <  N  )  =  1  . 

CONTINUE 
GO  TO  700 
X  -  4. L07 
MGF  <  N  >  =  t 
1  HI  (  N  >  / 


*  DIF**3) 


'  LAMBDA 
,16. 283 
AMBHA 


*  < 1  /  (DELTA  *  X#*C1 > 


l./X)  * 


HOF (N >  =  HGF < N )  -  39.48  *  HT(N>**2  *  (J.  X**C2  * 

1  HFL  FA  *  TAO)  /  <  PEL  TA  *  L  AMBDA**2  *  X**C3> 

CONT I NOE 

SUMX  a  Hlif  <  1  >  *  H(Jf  (2)  *  EXP<J  *  TAU  *  CHI) 
l  /  <2.  *  TAU  -  1.  /  DEL  1  A**2 ) 

SUMX  -  SUMX  *  fc 
SUM  *  SUM  f  SUMX 

IF  <  1  .  E  -  4  *  AHS(SOM) .GT. ASS (SUMX) >  GO  10  900 

con  r 1 NOt 

E  ABS(SUM) 

******  ****** ************************** 

COMMIT  fc  RECEIVED  POWER  IN  DDU 


************************************** 

IF  (LAMBDA.!  E.O. >  LAMBDA  =  1. 

IF  (E.LT.O. )  E  «  1 . 

IF  < PUR .11.0.)  PUR  =  1. 

RFCPWR  =  -158.513  f  20.  *  l  GIG  1 0 < LAMBDA )  f  20.  *  L0G10<E)  f 
1  10.  *  lOGlO(PWR) 

************************************** 


OUT  PUT  SECTION 


C  ************************************** 

SIGMA  =  SIGMA  *  l.EU 

UR 1 TE <  2  r 1 000  >  E r  RECPWR  >  PUR  r  LAMBDA » D t HT  » NPOL  r  FPSLN  *  SIGMA 
1000  FORMAT  </, 1X» 'OUTPUT: 7 .///» 1X» 'FIELD  STRENGTH  AT  R  ■  '»4X» 

1  IF’E'V. 3f  1X»  'MICROVOLTS/METER'  r//» IX*  'POWER  AT  R  5  '*7X*OPF10.2 
2 1 6X  f ' DBU ' ,///*  tX» '************************* #  t///t IX* ' INPUT  *  '  * 
3///»IX»'T  FOUFR  -  '*11X*F9.3»6X*' KILOWATTS' * 

4/  /  *  1 X  *  'WAVELENGTH  =  '  »8X»F9.3»6X»  'METERS'  * // *  IX  *  'DISTANCE 
5  MR  =  ' »/X/f0.3r6X, 'KILOMETERS' *//* IX* 'HEIGHT  OF  T  ■  '*7X 
6 » f-  / ,  1  » 8X »  ' ME  T ERS ' * // *  1 X *  '  ME  I GHT  OF  R  =  '*  8X  *F7 . 1  *  0X  •  '  METERS '  » 
7//*lX* 'POLARIZATION  *  ' 1 1 0X * 1 1  * T3B* ' 0  =  VERTICAL  *  1  = 

H  H0RIZ0NTA1  '»//* 1X»  EPSILON  =  '  ? 1 4X  »F  5  .2*//*lX*' SIGMA  =  '*17 
9X  *  1  PE 8 . 2 1 3X * 'MHOS/ME  IER' ) 

END 
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Figure  1 


Comparison  of  ground  wave  r.m.s.  field  strengths  calculated  by 
the  MRC  ground  wave  program  (x)  and  by  H.  Bremmer.1  The 
ground  constants  apply  to  sea  water,  the  polarization  is  verti¬ 
cal,  the  transmitter  power  is  1  kw,  and  both  transmitter  and 
receiver  are  on  the  ground. 
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Figure  3.  Comparison  of  basic  transmission  loss  in  dB  calculated  by  the 
MRC  ground  wave  program  (x)  and  by  K.  Norton.12  The  compu¬ 
tations  are  for  average  land  (o  =  .005  mhos/meter,  e  =  15), 
vertical  polarization  and  both  transmitting  and  receiving 
antennas  nine  meters  above  the  ground. 


accomplished  with  an  earlier  version  of  the  program  using  Bremmer's  origi¬ 
nal  conventions  for  transmitter  power  and  distance  between  transmitter  and 
receiver,  not  the  CCIR  conventions  now  used.)  Figures  1  and  2  show  that 
the  MRC  ground  wave  program  matches  Bremmer's  computations  almost  per¬ 
fectly.  The  very  small  differences  are  probably  due  to  the  greater  accur¬ 
acy  possible  with  a  computer.  (All  of  Bremmer's  results  were  computed  by 

hand.)  In  contrast,  the  comparisons  in  Figure  3  between  MRC  and 
Norton^  reveal  small  although  significant  discrepancies.  MRC's  results 
are  consistently  lower  than  Norton's.  These  discrepancies  probably  arise 
because  Norton  also  considers  atmospheric  refraction  while  we  do  not. 
Refraction  is  discussed  further  in  Section  4. 

The  MRC  results  also  match  the  world  standard  CCIR  results  from 
the  Geneva  conference  of  1974, but  are  somewhat  lower  than  the  CCIR 
results  from  the  Kyoto  conference  of  1978. 14  This  is  not  surprising 

because  the  1974  results  were  based  on  Bremmer's  equations  while  the  1978 

results  were  based  on  Bremmer's  equations  but  also  included  refraction.  A 
set  of  CCIR  ground  wave  propagation  curves  is  provided  for  reference  in 
the  appendix. 

The  approach  used  by  the  MRC  program  is  quite  similar  to  that 
used  by  the  NUCOM/BREM  program.^ 

Using  computations  from  the  MRC  ground  wave  program,  the  curves 
in  Figures  4,  5,  and  6  were  generated.  These  curves  apply  to  a  1  kW 
transmitter  and  to  wide  variety  of  circumstances  from  sea  water  to  average 
land,  from  a  frequency  of  15  KHz  to  10  GHz,  from  a  range  of  2  km  to  2000 
km,  and  from  a  receiver  height  of  0  m  to  1000  m.  They  demonstrate  the 
versatility  of  the  program  while  not  attempting  to  cover  all  possible 
parametric  variations.  Figure  6  also  shows  that  the  program  fails  to 
compute  correct  values  for  a  high  receiver  height  and  short  distance 
between  transmitter  and  receiver.  The  departures  from  a  smooth  curve  at 
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Figure  5. 


-40 


Figure  6.  Similar  to  Figures  4  and  5  except  the  frequency  here  is  fixed 
at  30  MHz  and  the  height  of  the  receiver  (in  m)  is  varied.  The 
transmitter  height  is  fixed  at  10  m  and  the  curves  apply  to 
average  ground. 


receiver  heights  of  400  and  1000  m  are  errors.  The  program  generally 
fails  at  very  short  separations  between  transmitter  and  receiver.  In 
these  cases,  a  1 ine-of-sight  or  geometric  optical  model  is  more  appro¬ 
priate. 


Figures  7  and  8  show  one  application  of  the  MRC  ground  wave 
program  to  a  hypothetical  ground  wave  system.  The  ground  constants  used 
are  appropriate  to  average  ground  in  Germany  while  the  transmitter  charac¬ 
teristics  used  are  typical  of  many  commercially  available  products.  The 
maximim  range  is  determined  by  the  noise  level  at  the  receiver  and  not  the 
receiver  itself.  Two  curves  are  shown  in  each  figure  corresponding  to  a 
best  and  worst  case  noise.  The  actual  situation  will  probably  be  some¬ 
where  in  between  the  two  extremes.  Figure  8  shows  the  considerable 
improvement  in  range  that  is  possible  with  nominal  improvements  in  trans¬ 
mitter  characteristics. 
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Figure  7.  Maximum  ground  wave  range  for  a  hypothetical  system  In  Germany. 
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Figure  8.  Same  as  Figure  7  except  for  improved  transmitter  characteris 
tics. 
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SECTION  4 
REFRACTION 


This  section  discusses  the  refraction  of  ground  waves  by  a 
non-uniform  atmosphere.  Figure  9  shows  that  refraction  is  a  significant 
effect  which,  under  most  but  not  all  circumstances,  improves  ground  wave 
propagation. 

The  variation  of  refractive  index  as  a  function  of  altitude  can 
be  expressed  in  several  ways.  The  following  equations  show  two 
possibilities 


and 


p(h)  =  1  +  A  exp  (-Bh) 


(19) 


(20) 


where  A  and  B  as  well  as  n  and  y  are  arbitrary  constants  chosen  to  suit 
the  particular  location  involved,  h  is  the  height  above  mean  sea  level, 
r  is  the  distance  from  the  center  of  the  earth,  and  a  is  the  radius  of  the 
earth.  Equation  19  comes  from  the  CCIR*4  while  Equation  20  comes  from 
Bremmer.l  Although  19  is  now  considered  the  preferable  form,  20  is 
adequate  in  the  lower  atmosphere.  We  choose  to  use  Equation  20  because 
Bremmer's  equations  for  refraction  readily  follow  from  this  equation. 

Introducing  the  parameter  a. 


a 


1-n 

1-n+y 


(21) 
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it  can  be  shown*  that  the  electric  field  E  at  a  distance  D  and  height  h2 
from  a  transmitter  at  height  ht  over  soil  with  a  parameter  6  (see  Equation 
9)  can  be  reduced  to  the  electric  field,  E0,  that  would  exist  without 
refraction  in  the  following  way: 

E(D,h1,h2,6)  =  ot2/3  Eo(Da2/3,  h^1/3,  h261/ 3,  6a1/3)  (22) 

In  other  words.  Equation  6  can  be  used  to  compute  the  electric  field  with 
refraction  providing  the  modified  parameters  shown  in  Equation  22  are 
used.  It  is  also  necessary  to  use  the  actual  wavelength  existing  at  the 
earth's  surface  x0/u(a)  and  not  the  wavelength  in  vacuum. 

To  be  able  to  use  Equation  22,  it  is  necessary  to  evaluate  a. 
This  can  be  done  using  the  following  formula  for  the  index  of 
refraction:*!*4 


where 


u  =  1  +  77.6  E.  +  .373  -1 
T  T2 

T  =  absolute  temperature  (’K) 
p  =  atmospheric  pressure  in  millibars 
e  =  water  vapor  pressure  in  millibars 


Differentiating  20  and  23  with  respect  to  height  yields. 


and 


u'(r)  =  - 


T 

au(r) 


(23) 


(24) 


ii'(h)  =  -77.6*10-6  _P  £L  +  7.7--.6.x10.'. .  jR.  -  .746  —  —  +  : i73.  —  (25) 

T2  dh  T  dh  T3  dh  T2  dh 


Equation  25  can  be  simplified  by  use  of  the  relation 
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(26) 


dp  de 
dh  dh 


-10  gp 


where  the  unit  of  length  is  chosen  to  be  100  m  and 

g  =  acceleration  of  gravity  in  cm/sec2 
p  =  density  of  air  in  g/cm2. 

Substituting  Equation  26  in  25  yields 

y '  (h)  =  -  jZZ  p  +  f-77.6xl0-6  +  —1 ) 

T  T  T  dh 


(27) 


-  (77.6xl0“6  p  +  .746  -)  —  — 

T  T2  dh 

Because  refraction  takes  place  close  to  the  earth's  surface,  p(r)  and 
p'(r)  can  be  adequately  approximated  by  conditions  at  the  earth's 
surface.  Because 

Y-n  =  u2(a)-l  ■  .00058  *  0  (28) 

Y-n  can  be  neglected.  Then  (24)  can  be  written 

u'(a)  =  -  —  (29) 

a 

Equations  21,  27,  and  29  can  be  combined  to  provide  the  following  solution 
for  a  in  terms  of  atmospheric  conditions  near  the  earth's  surface 

a  -  1  -  a[-lZ®  p  +  (77.6xl0-6  +  (30) 


+  (77.6xl0_6  p  +  .746  -)  —  — 

T  T2  dh 
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Equation  30  when  used  in  conjunction  with  Equation  22  is  an  adequate 
formal  solution  of  the  refraction  problem  which  can  be  readily  adapted  for 
use  on  a  computer.  Using  a  better  approximation  than  Equation  29  will 
lead  to  a  slightly  better  solution  of  the  refraction  problem. 

Without  a  computer,  the  formal  solution  of  Equations  22  and  30 
can  still  be  used  by  employing  several  additional  approximations.  Using 
the  approximation  in  Equation  28,  Equation  22  can  be  rewritten 

E(D,h x ,h2)  =  (1-n)4/3  Eo[n(l-n)2/3,  h^l-n)1/3,  h2(l-n)1^3]  (31) 

This  equation  further  assumes  that  6  (given  in  Equation  9)  is  small.  That 
is, 

| 6f  =  K  «  1  .  (32) 

This  occurs  for  short  wavelengths  from  VHF  and  shorter  (perhaps  also  for 
portions  of  the  HF  band  depending  on  the  ground  constants  used).  In  any 
case.  Equation  31  is  considerably  simpler  than  22. 

Without  refraction  the  electric  field  depends  on  D  principly 
thru  the  parameter,  given  in  Equation  7: 

X  =  (2ir)l/3  D/(az/3  Xl/  3)  =  .053693  D^/^1/3  (33) 

From  this  expression,  it  is  evident  that  changing  D  to  D(l-n)2^3  is 
equivalent  to  changing  the  radius  of  the  earth  to  a/(l-n)  and  holding  D 

constant.  The  effective  radius  of  the  earth  can  then  be  written 

a  =  a  /  [  .766  -  (.0681  +  .00237  e)  —  +  .306  —]  (34) 

e  L  dh  dh J 
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using  Equations  27  and  29  and  evaluating  T  and  p  at  standard  temperature 
and  pressure  (STP). 

Equation  34  is  a  most  interesting  expression.  It  is  evident  that 
for  a  dry  atmosphere  with  an  adiabatic  temperature  gradient  (dT/dh  «-l), 
ae  «  1.20a.  For  saturated  air,  ae  «  1.55a.  For  mean  meteorological 
conditions  ae  *>  (4/3)a,  which  is  an  often  used  value.  Under  some  condi¬ 
tions  ae  <  a,  and  the  refraction  will  actually  hinder  ground  wave  propa¬ 
gation.  In  contrast,  a  temperature  inversion  (dt/dh  >  0)  will  considerably 
enhance  radio  propagation. 

To  compute  the  resulting  electric  field  in  the  case  of  refrac¬ 
tion  using  the  approximations  above,  it  is  only  necessary  to  take  a  graph 
of  E  versus  D  without  refraction  and  multiply  E  by  (1-n)*4/3  and  D  by 
l/(l-n)2/3.  If  hL  and  h2  are  non  zero,  then  they  must  be  similarly 
reduced  as  shown  in  Equation  31.  The  graph  in  Figure  9  was  prepared  by 
Reference  1  using  this  technique. 


SECTION  5 

PROPAGATION  OVER  IRREGULAR  TERRAIN 


Up  to  this  point,  we  have  assumed  that  the  earth  is  a  smooth 
homogeneous  sphere.  Although  the  earth  varies  substantially  from  this 
idealization,  the  results  expressed  in  Equation  6  are  applicable  to  many 
situations.  Nevertheless,  it  would  be  useful  to  compute  field  strengths 
for  those  cases  where  Equation  6  is  clearly  not  valid.  For  instance, 
propagation  over  mountains  or  propagation  from  land  to  ocean  is  not  well 
described  by  Equation  6.  The  literature  on  these  subjects  is  extensive 
and  well  beyond  the  scope  of  this  paper.  We  arbitrarily  choose  to  discuss 
propagation  over  obstacles  and  to  refer  the  reader  to  other  work  for 
discussions  of  propagation  over  inhomogeneous  or  inductive  terrain.  In 
all  of  these  cases,  the  results  occasionally  differ  from  simple  intuitive 
extensions  of  Equation  6  because  of  the  unusual  nature  of  diffraction 
processes.  For  instance,  a  large  obstacle  between  a  transmitter  and 
receiver  can  actually  increase  the  received  signal  strength.  Such  an 
effect  is  called  "obstacle  gain"  and  is  explained  below. 

Propagation  over  inhomogeneous  terrain  can  be  handled  in  various 
ways.  The  Suda  method^  involves  averaging  the  ground  constants  for  var¬ 
ious  segments  of  the  terrain  while  the  Millington  method^7*!8  involves  a 
geometrical  mean  of  the  electric  fields  due  to  the  various  segments.  The 
Suda  method  is  most  appropriate  where  the  ground* constants  do  not  change 
substantially  while  the  Millington  method  is  applicable  to  substantial 
changes  as  long  as  the  field  is  measured  well  away  from  the  interface 
region.  Many  other  references  have  considered  these  approaches  as  well  as 
others.  References  15  and  19  through  31  discuss  the  subject  in  detail. 
In  addition.  Reference  32  discusses  propagation  over  inductive  terrain. 
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Propagation  over  obstacles  has  also  been  covered  extensively  in 
the  literature.  The  problem  can  be  expressed  in  a  formal  manner  in  terms 
of  a  two  dimensional  integral  equation  of  the  following  type:^»^3 


W(p) 


1 


J. 

2tt 


//  W(Q)  eiko^TQ  +  QP  -  PT) 

QPxQT 


x  [y(Q)  -  (ik  -  L)W]  dO 
QP  anQ 


(35) 


U  is  the  ratio  of  the  actual  amplitude  of  the  Hertzian  vector  to  its  value 
for  the  case  of  a  perfectly  conducting  plane  earth.  The  distances  TQ,  QP, 
etc.  are  indicated  in  Figure  10.  The  function  y  is  related  to  the  tilt  of 
the'  field  as  well  as  to  so-called  surface  admittances.  The  integral 
extends  over  the  entire  irregular  portion  of  the  earth's  surface  except 
for  an  infinitesimal  area  around  P.  This  equation  can  account  for  both 
irregular  and  inhomogeneous  terrain. 


Equation  35  can  be  reduced  to  a  one-dimensional  equation 
provided  that  the  irregular  surface  of  the  earth  does  not  deviate  too 
greatly  from  an  average  level  plane  or  that  the  various  ground  parameters 
do  not  vary  too  greatly  perpendicular  to  the  line  connecting  the 
transmitter,  T,  and  receiver,  P.  A  saddlepoint  approximation  simplifies 
35  to 


Wt(x) 


jx  -MlL 

0  /C(x- 5) 


[yiU)  +  ik0 


sin 


«l(0] 


eik0(TQ+QP-PT)dc 


(36) 


This  expression  shows  that: 

a)  Local  changes  of  the  ground  constants  and  of  the  terrain 
profile  have  a  similar  effect  on  Wj.  The  term  in  square 
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Figure  10.  Geometrical  representation  of  an  irregular  terrain  profile. 
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brackets  shows  that  y»  representing  ground  constants, 
and  a,  representing  terrain  profile,  have  a  similar 
effect. 


b)  The  terrain  near  the  transmitter  and  receiver  has  a 
relatively  greater  influence  due  to  the  existence  of  the 
weighting  factor  l//£(x-c). 

It  is  possible  to  reduce  Equation  35  to  Equation  36  because  of 
the  dominating  role  of  integration  points,  Q,  near  stationary  values  of  an 
exponential  phase  factor.  The  region  surrounding  the  area  of  a  stationary 
phase  is  the  first  Fresnel  zone.  Figure  11  shows  a  drawing  of  this  zone. 
Terrain  irregularities  and  inhomogeneities  produce  significantly  different 
results  if  they  occur  inside  or  outside  this  first  Fresnel  zone.  Inside 
the  zone,  disturbances  can  prevent  production  of  the  main  field  at  the 
receiver.  This  is  the  field  that  would  have  occurred  had  no  disturbances 
been  present.  In  contrast,  disturbances  outside  the  zone  will  not  signif¬ 
icantly  alter  the  main  field  but  may  add  some  perturbations.  If  the 
transmitter  and/or  receiver  are  elevated,  the  same  reasoning  applies  but 
the  first  Fresnel  zone  is  somewhat  different. 

Major  obstacles  between  the  transmitter  and  receiver,  such  as 
the  ridge  shown  in  Figure  12,  can  produce  a  surprising  result  known  as 
"obstacle  gain".  The  ridge  is  assumed  to  be  within  the  first  Fresnel  zone 
and  well  above  the  shortest  path  between  the  transmitter  and  receiver.  In 
this  case,  the  ridge  acts  approximately  like  the  absorbing  straight  edge 
often  considered  in  the  theory  of  optical  diffraction. 2  The  field  at  or 
beyond  the  ridge  can  be  estimated  from  the  Cornu  spiral  shown  in  Figure 
13.  The  starting  point  for  the  spiral  is  fixed  at  the  lower  point  of  the 
spiral,  and  the  various  chords,  a  through  g,  correspond  to  the  magnK'de 
and  phase  of  the  electric  field.  Only  the  end  point  of  the  chord  changes 
with  x.  (x  is  the  coordinate  perpendicular  to  the  straight  edge  and 
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parallel  to  the  wave  front.)  In  the  geometrical  shadow  region  (-»  <  x  <  0) 
the  length  of  the  chord  increases  steadily  as  indicated  by  points  a 
through  d.  d  corresponds  to  the  boundary  of  the  shadow  region.  From  this 
point,  the  chord  keeps  increasing  in  length  until  it  reaches  a  maximim  at 
f  and  then  begins  an  oscillatory  behaviour  approaching  the  other  end  of 
the  sprial.  Figure  14  shows  this  effect  explicitly. 

Behind  the  diffracting  straight  edge  the  field  is  given  approxi¬ 
mately  by 


-  + -  U/; 

cos(l^)  cos(£H) 

where  u  represents  the  electric  field  when  it  is  parallel  to  the  edge  and 
the  minus  sign  applies.  When  the  electric  field  is  perpendicular  to  the 
edge,  u  represents  the  magnetic  field  and  the  plus  sign  applies.  Equation 
37  assumes  a  cylindrical  coordinate  system  with  the  z-axis  parallel  to  the 
straight  edge  and  the  angular  coordinate  <f>=0  or  2n  at  the  straight  edge, 
a  is  the  angle  between  the  wave  normal  and  straight  edge.  Because  the 
value  in  square  brackets  decreases  slowly  with  increasing  «)>,  the  light  or 
radio  wave  is  diffracted  far  into  the  shadow  region.  Equation  37  is  based 
on  an  approximation  which  fails  at  the  shadow  boundary  <j>  =  a  or  $  =  -a. 

Equation  37  shows  that  the  diffracted  field  beyond  the  ridge 
falls  off  as  r-1/2.  This  is  less  of  a  decrease  than  for  ground  wave 
propagation  over  a  flat  earth  (r-2)  or  a  curved  earth  [e~r^r°).  In 
front  of  the  ridge,  the  field  decreases  as  r_1  which  is  also  better  than 
ground  wave  propagation.  In  other  words,  a  sharp  ridge  can  actually 

improve  reception  at  a  distant  location  behind  the  ridge.  Because  typical 
soil  conditions  produce  considerable  absorption  of  ground  waves  while  a 
propagation  path  over  a  high  ridge  avoids  the  ground  to  a  large  extent, 
improved  propagation  is  possible  by  utilizing  a  ridge. 


(1+i)  gikr 


4/ttTct' 
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Figure  13. 


Figure  14. 


Cornu  spiral  construction  for  determining  the  diffraction 
pattern  produced  by  a  straight  edge.  The  axes  labeled  S  and  C 
represent  the  Fresnel  integrals. 


Amplitude  of  the  electric  field  at  a  straight  edge.  Zero 
corresponds  to  the  boundary  of  the  edge*  negative  values  of  w 
are  behind  the  edge  Mhile  positive  are  in  front  of  the  edge. 
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Several  computer  programs  have  been  developed  to  compute  ground 
wave  propagation  over  irregular  and  inhomogeneous  terrain.  These  require 
a  detailed  knowledge  of  the  terrain  between  the  transmitter  and  receiver 
and  involve  long  or  short  wavelength  approximations.  One  program  devel¬ 
oped  by  the  Communications  Research  Centre-^4  applies  to  frequencies  at 
VHF  or  higher.  Another,  developed  by  ITS^S  applies  to  MF  and  lower  fre¬ 
quencies  but  can  be  stretched  into  the  HF.  The  basic  difficulty  with  the 
HF  is  that  the  radio  wavelength  is  close  to  the  size  of  terrain  features 
so  that  both  long  and  short  wavelength  approximations  have  limited  applic- 
abi 1 ity. 


Another  difficulty  in  the  HF  band  is  the  decided  lack  of  experi¬ 
mental  data  to  help  verify  computer  routines.  While  ground  wave  data  for 
rough  terrain  exists  at  MF  and  VHF,  it  is  surprising  that  little  exists  at 
HF.  Perhaps  the  HF  has  been  ignored  because  it  is  used  principly  for 
longer  distance  sky  wave  propagation.  In  any  case,  data  at  HF  is  essen¬ 
tial  to  further  our  ability  to  accurately  predict  HF  propagation  over 
rough  terrain. 

Several  investigators  have  collected  data  over  rough  terrain  at 
frequencies  outside  the  HF.^-4^  gf  these  the  work  by  ITS  is  perhaps 
the  most  thorough. 35,36  xhis  da-t a  has  teen  used  to  show  that  the  ITS 
program  called  WAGNER  is  indeed  able  to  predict  ground  wave  propagation 
over  rough  terrain  below  a  few  megahertz. 
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APPENDIX 


The  following  curves,  considered  to  be  the  best  available 
predictions  of  ground  wave  field  strengths,  are  reproduced  from  the  CCIR 
conference  of  1978  in  Kyoto,  Japan.  They  come  from  the  corrigendum  to 
volume  V,  Propagation  in  Non-Ionized  Media,  and  are  dated  February  25, 
1980. 


The  curves  refer  to  the  following  conditions: 


a)  smooth  homogeneous  earth 

b)  transmitter  and  receiver  are  on  the  earth 

c)  transmitting  antenna  is  an  ideal  Hertzian  vertical 
electric  dipole  (nearly  identical  to  a  vertical  antenna 
shorter  than  one  quarter  wavelength) 

d)  transmitter  power  is  defined  in  exactly  the  same  way  as 
indicated  in  Section  3  (basically  1  kw). 

e)  refraction  is  accounted  for  by  assuming  an  atmosphere  in 
which  the  refractive  index  decreases  exponentially  with 
height. 

f)  ionospheric  reflections  are  excluded. 

g)  distances  are  measured  around  the  curved  surface  of  the 
earth. 
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Inverse  distance  curve 


Ground  wave  propagation  curves;  Land 


Inverse  distance  curve 
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ATTN:  J.  Hirsch 

Boeing  Co 

ATTN:  S.  Tashird 
ATTN:  G.  Hall 

ATTN:  M/S  42-33,  J.  Kennedy 

Booz-Allen  A  Hamilton,  Inc 
ATTN:  B.  Wilkinson 

University  of  California  at  San  Diego 
ATTN:  H.  Booker 

Charles  Stark  Draper  Lab,  Inc 
ATTN:  J.  Gilmore 
ATTN:  D.  Cox 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Computer  Sciences  Corp 

ATTN:  F.  Eisenbarth 

Comsat  Labs 

ATTN:  D.  Fang 
ATTN:  G.  Hyde 

Cornell  University 
ATTN:  M.  Kelly 
ATTN:  D.  Farley,  Jr 

E-Systems,  Inc 

ATTN:  R.  Berezdivin 

Electrospace  Systems,  Inc 
ATTN:  H.  Logston 

ESL,  Inc 

ATTN:  J.  Lehman 
ATTN:  R.  Ibaraki 
ATTN:  R.  Heckman 
ATTN:  J.  Marshall 
ATTN:  E.  Tsui 

General  Electric  Co 

ATTN:  G.  Holliday 
ATTN:  A,  Harcar 

General  Electric  Co 

ATTN:  A.  Stei nmayer 
ATTN:  C.  Zierdt 

General  Electric  Co 

ATTN:  F.  Reibert 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
JAYCOR 


ATTN 

J.  Sperling 

JAYCOR 

ATTN 

J.  Doncarlos 

Johns 

Hopkins  University 

ATTN 

T.  Potemra 

ATTN 

J.  Phillips 

ATTN 

T.  Evans 

ATTN 

J.  Newland 

ATTN 

P.  Komiske 

Kaman 

Sciences  Corp 

ATTN 

T.  Stephens 

Kaman 

Tempo 

ATTN 

DAS  I AC 

ATTN 

W.  Knapp 

ATTN 

W.  McNamara 

Linkabit  Corp 

ATTN 

I.  Jacobs 

ATTN 

A.  Viterbi 

ATTN 

H.  Van  Trees 

Litton  Systems,  Inc 

ATTN 

R.  Grasty 

Lockheed  Missiles  &  Space 

ATTN 

W.  Imhof 

ATTN 

R.  Johnson 

ATTN 

M.  Walt 

General  Electric  Co 


ATTN: 

G. 

Mil lman 

General  Research  Corp 

ATTN: 

J. 

Garbarino 

ATTN: 

J. 

Ise,  Jr 

Harris  Corp 

ATTN: 

E. 

Knick 

Horizons  Technology,  Inc 
ATTN:  R.  Kruger 

HSS,  Inc 

ATTN:  D.  Hansen 
IBM  Corp 

ATTN:  H.  Ulander 

University  of  Illinois 
ATTN:  K.  Yeh 


Institute  for  Defense  Analyses 


ATTN: 

H.  Gates 

ATTN: 

J.  Aein 

ATTN: 

E.  Bauer 

ATTN: 

H.  Wolfhard 

Internationa 

Tel  &  Telegraph  Corp 

ATTN: 

W.  Rice 

ATTN: 

Technical  Library 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  C.  Old,  Dept  68-21 
ATTN:  D.  Churchill,  Dept  81-11 
ATTN:  Dept  60-12 

M.I.T.  Lincoln  Lab 

ATTN:  D.  Towle 

Magnavox  Govt  &  Indus  Electronics  Co 
ATTN:  G.  White 

Martin  Marietta  Corp 

ATTN:  R.  Heffner 

McDonnell  Douglas  Corp 
ATTN:  H.  Spitzer 
ATTN:  W.  Olson 
ATTN:  R.  Halprin 

Meteor  Communicatons  Consultants 
ATTN:  R.  Leader 


Mission  Research  Corp 


ATTN 

R.  Kilb 

ATTN 

Tech  Library 

ATTN 

D.  Sappenfield 

ATTN 

F.  Fajen 

ATTN 

R.  Hendrick 

ATTN 

R.  Bogusch 

ATTN 

S.  Gutsche 

4  cy  ATTN 

G.  Fulks 

5  cy  ATTN 

Doc  Con 
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Mitre  Corp 

ATTN 

B.  Adams 

ATTN 

G.  Harding 

ATTN 

A.  Kymmel 

ATTN 

C.  Callahan 

Mitre  Corp 

ATTN 

M.  Horrocks 

ATTN 

W.  Foster 

ATTN 

W.  Hall 

ATTN 

J.  Wheeler 

Pacific-Sierra  Research  Corp 
ATTN:  F.  Thomas 
ATTN:  E.  Field,  Jr 
ATTN:  H.  Brode,  Chairman  SAGE 

Pennsylvania  State  University 

ATTN:  Ionospheric  Research  Lab 

Photometries,  Inc 

ATTN:  I.  Kofsky 

Physical  Dynamics,  Inc 
ATTN:  E.  Fremcuw 

Physical  Research,  Inc 


ATTN 

R.  Deliberis 

R  &  D  Associates 

ATTN 

F.  Gilmore 

ATTN 

B.  Gabbard 

ATTN 

M.  Gantsweg 

ATTN 

W.  Karzas 

ATTN 

W.  Wright 

ATTN 

C.  Greifinger 

ATTN 

R.  Lelevier 

ATTN 

H.  Ory 

A1TN 

R.  Turco 

ATTN 

P.  Haas 

R  &  D  Associates 

ATTN 

:  B.  Yoon 

Rand  Corp 

ATTN:  E.  Bedrozian 
ATTN:  C.  Crain 

Riverside  Research  Institute 
ATTN:  V.  Trapani 

Rockwell  International  Corp 
ATTN:  R.  Buckner 


Santa  Fe  Corp 

ATTN:  D.  Paolucci 


Science  Applications,  Inc 

ATTN 

C.  Smith 

ATTN 

L.  Linson 

ATTN 

D.  Hamlin 

ATTN 

E.  Straker 

Science  Applications,  Inc 

ATTN 

SZ 

Science  Applications,  Inc 

ATTN 

J.  Cockayne 

SRI  International 

ATTN 

A.  Burns 

ATTN 

R.  Livingston 

ATTN 

R.  Leada brand 

ATTN 

D.  Neilson 

ATTN 

M.  Baron 

ATTN 

W.  Chesnut 

ATTN 

R.  Tsunoda 

ATTN 

G.  Smith 

ATTN 

G.  Price 

ATTN 

C.  Rino 

ATTN 

W.  Jaye 

ATTN 

J.  Petrickes 

Sylvania  Systems  Group 

ATTN 

J.  Concordia 

ATTN 

I.  Kohl  berg 

ATTN 

R.  Steinhoff 

Technology  International  Corp 

ATTN:  Doc  Con  for  W.  Boquist 

Tri-Corn,  Inc 

ATTN:  D.  Murray 

TRW  Electronics  &  Defense  Sector 
ATTN:  R.  Pleubuch 
ATTN:  D.  Dee 

Utah  State  University 
ATTN:  K.  Baker 
ATTN:  L.  Jensen 
ATTN:  J.  Dupnik 

Visidyne,  Inc 

ATTN:  C.  Humphrey 
ATTN:  J.  Carpenter 


Rockwell  International  Corp 
ATTN:  S.  Quilici 


